Distribution of Interference in the Presence of Decoherence 



Ludovic Arnaud and Daniel Braun 

Universite de Toulouse, UPS, Laboratoire de Physique 
Theorique (IRSAMC), F-31062 Toulouse, France and 
CNRS, LPT (IRSAMC), F-31062 Toulouse, France 

Abstract 

We study the statistics of quantum interference for completely positive maps. We calculate 
analytically the mean interference and its second moment for finite dimensional quantum systems 
interacting with a simple environment consisting of one or several spins (qudits). The joint propa- 
gation of the entire system is taken as unitary with an evolution operator drawn from the Circular 
Unitary Ensemble (CUE). We show that the mean interference decays with a power law as func- 
tion of the dimension of the Hilbert space of the environment, with a power that depends on the 
temperature of the environment. 
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I. INTRODUCTION 



Quantum information theory predicts increased computational power for quantum algo- 
rithms compared to classical algorithms. The most well-known example is Shor's algorithm 
which factors a large integer number in a time which grows only polynomially in the num- 
ber of digits [1], whereas no such algorithm is known classically. Grover found a quantum 
algorithm that allows to find an item in an unstructured data base of size N with a number 
of queries that scales only like VN, whereas classically the number of queries is of order N. 
Exponential acceleration compared to the best known classical algorithm was also predicted 
for the shifted character problem the hidden subgroup problem 3], and for solving linear 
systems of equations 4j. A quantum walk can traverse a graph exponentially faster than 
any classical random walk which allows for the efficient solution of certain oracle problems 
j^]. Aharonov et al. proposed a quantum algorithm which efficiently approximates the Jones 
polynomial at any primitive root of unity |6j. 

It seems to be clear that quantum entanglement and quantum interference are two key 
resources which provide for the enhanced information processing capabilities of quantum 
systems Q. But in spite of the many known examples in which quantum information 
processing outperforms classical information processing, it is not entirely clear how exactly 
these resources enable the speed of quantum algorithms, nor what the largest possible speed- 
up is. It was shown [8| that a unitary quantum algorithm in which entanglement remains 
"p-blocked" (i.e. the number of qubits which at any time are entangled is not larger than p), 
can be efficiently simulated classically. Nevertheless, the same authors argued that it might 
be misleading to consider entanglement as the key resource. As long as the mechanism is 
not identified by which any specific quantity creates the speed-up, one might suspect its 
creation in large amounts rather correlated with the quantum acceleration than being its 

CI 

cause. Entanglement is definitely crucial for tasks like quantum teleportation [9j, where its 
role can be understood through the enhanced correlations between subsystems that quantum 
mechanics can provide. 

Recently, experimental implementations lOj of factoring integer numbers using Gauss 
sums, have re-emphasized the role of interference in quantum computation. While these 
methods do not appear to be scalable to integers with many 100 digits, and can be im- 
plemented with classical waves, they are reminiscent of simple quantum algorithms like the 
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Deutsch-Jozsa algorithm, in which interference is clearly seen at work. Contrary to quantum 
entanglement, quantum interference has been surprisingly little studied. From a physicist's 
perspective, quantum interference is an effect that arises from the coherent superposition 
of quantum mechanical wave functions. This can lead to interference maxima and min- 
ima in probability distributions, as is well-known from quantum particles going through a 



double-slit, electrons in a mesocscopic solid state circuit [ll||, or interfering Bose-Einstein 
condensates [12| . Quantum interference can also focus the probability distribution in a com- 
puter over its possible states at the outcome of a calculation onto the state corresponding 
to the result of the calculation. Without the coherence of quantum superpositions, proba- 
bilities can only be propagated classically, i.e. through a stochastic map, which is, of course, 
void of any interference effects. If we want to quantify interference, we therefore have to 
quantify to what extent the propagation is coherent, as otherwise there is no telling if the 
production of a final probability distribution involved interference or not. But coherent 
propagation alone is not tantamount to interference. At least two wave functions have to be 
superposed in order to create interference. Very generally, one would want to attribute more 
interference to a process in which many waves get superposed with similar weights than to 
one where only very few waves contribute. This implies a basis dependence of interference, 
as a superposition in one basis is a single basis state in another. 



In 13j a measure of quantum interference was introduced which allows to quantify inter- 
ference in any quantum mechanical process in a finite dimensional Hilbert space. Any such 
process can be described by a completely positive map P that maps an initial density matrix 
p to a final one, p' = Pp. Written in a given basis, where p and p' have matrix elements p mn 
and p' kl , respectively, we have p' mn = ^2 kl P mn kiPki- In that basis, the interference associated 
with the positive map P is written as 

%{P) = S ^\Pii,kl\ 2 ~' S ^\Pii ) kk\ 2 ■ (1) 

While this interference measure may not be unique, it has the desired property of measur- 
ing the coherence and the "equipartition" of superposed basis states. Indeed, if P reduces 
to a classical stochastic map, it only propagates initial probabilities pa to final ones, p^k- 
Exactly the terms responsible for this classical process are subtracted out in eq.flTJ, such 
that if no coherences are propagated to final probabilities, we have zero interference. The 
squares of the matrix elements of P in ([I]) allow to measure the equipartition property, as 



is seen most easily for purely unitary propagation, where I reduces to iV — £V |£4/| 4 , where 
U is the unitary matrix propagating the wave function, and iV the dimension of Hilbert 
space. Perfectly equipartitioned unitary matrices (|Z7y| = 1/VN) create the maximum 
amount of interference possible for unitary propagation, X = N — 1. As an example, the 
Hadamard gate creates one bit of interference, an "i-bit". Both Shor's and Grover's algo- 
rithm create an exponential amount of interference (in the number of qubits). The part of 
the quantum algorithm after application of the initial Hadamard gates creates only about 
three i-bits in Grover's algorithm, but still an exponentially large amount of interference in 
Shor's algorithm. If the success probability of these algorithms is lowered by introducing 
unitary errors or decoherence, so is in general the interference 14j . For unitary quantum 
algorithms randomly drawn from the Circular Unitary Ensemble (CUE), interference is very 
narrowly distributed about the mean value, which itself is almost the maximum possible 



value 



15). In other words, almost all unitary quantum algorithms lead to an exponentially 



large amount of interference. This situation is reminiscent of entanglement, as almost all 
states of high-dimensional bipartite systems are close to maximally entangled 161 ] . 

It also turned also, however, that quantum interference is not necessary for several tasks. 
For example the transmission of a quantum state through a chain of qubits needs only a 
very small amount of interference [171 ] . And cloning of a quantum state can be performed 
just as well without interference as with interference fl8| . 

Almost all investigations of interference have focused so far on unitary propagation. Re- 
cently the benefits of more general, partly dissipative and decoherent evolutions have been 
emphasized, both in the context of quantum enhanced measurements 19], as in quantum 



computing 



201 ]. Moving on in this direction, we investigate in this paper the statistical 



properties of interference for general positive maps. We construct such maps by propa- 
gating unitarily a central system and an environment, which we take here both as finite 
dimensional quantum systems, and then tracing out the environment. We calculate analyt- 
ically the first and second moments of the distribution. We first focus on an environment 
that consists of a single spin (such ancillary qubit or qudit), and generalize then to 

an arbitrary number of spins, all taken initially in a thermal state at arbitrary temperature. 
We also calculate numerically the entire interference distribution for small system sizes. 
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II. STATISTICS OF INTERFERENCE FOR A QUANTUM SYSTEM COUPLED 
TO A SINGLE SPIN 



In this section we first review the propagation of a finite dimensional quantum system that 
interacts with an arbitrary environment consisting of another finite dimensional quantum 
system. The corresponding propagator is a completely positive map of the initial density 
matrix of the system to its final density matrix 2l| . While a finite dimensional environment 
does not constitute a true heat-bath in the sense of inducing irreversible behavior, the study 
of such a simple situation is motivated by quantum information theory, where one frequently 
encounters ancilla qubits that are added to the main quantum information processor. Fur- 
thermore, the tracing out of any environment with dimension larger than one does lead 
to decoherence as soon as the system and its environment become correlated or entangled, 
such that we will be able to study quantitatively the influence of decoherence on quantum 
interference. Further freedom lies in the choice of the initial state of the environment, which 
can be in a mixed state, e.g. a thermal state reached by interaction with its own heat-bath. 
We then derive the expression for the interference of a quantum system whose environment 
is a simple spin initially in thermal equilibrium, and study the statistical properties of the 
interference of the completely positive map of the system under joint unitary evolution of 
system and environment. 



A. Propagator for a completely positive map 

Consider a bipartite system consisting of a system S (Hilbert space Tis with dimension 
n) and an environment £ (Hilbert space "He with dimension m). Let W and W be the 
initial and final density matrices of the total system, respectively. We consider an initial 
product state W = a ® e of the density matrices a and e of the system and its environment, 
respectively. Under the condition that the total system "<S + £" can be considered closed 
on the time scale of the evolution we are interested in, the evolution of the system and its 
environment in the tensor product Hilbert space H.s^'He of dimension N = n x m is purely 
unitary and can be represented by a unitary matrix U, W = UWU^ . In components we 
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have 

n,m 

^a 1 a 2 ,/3 1 (3 2 = ^0:102,7172 ^71 llJlfoU Pifaji&z 

71)72,51,52 

where the indices with subscripts 1 and 2 label the basis states of the system and the 
environment, respectively. The final reduced density matrix of the system is found by 
tracing out the environment, p' = tr e W, or, explicitly, p' ai/3i = J2™ 2 W L ia 2,/3 ia2 ■ From (ED 
and the initial W 71J2 = p 7l< 5 1 e 72< 5 2 we obtain the propagation of S alone, 

n 

Pai/3i = ^ ] -Pqi/3i, 7 i5i P7i<5i 
7l>5l 

where the components of the propagator are given by 

fit 

Paf3,~f8 — U a ^^ v e up Uf i ^ &p (2) 

This propagator P is a superoperator that maps the initial density operator p to the final 
density operator p' . The procedure of "hamiltonian embedding" we have used guarantees 
that this propagator is a completely positive map [21]. As expected, P depends not only on 
U but also on the initial state of the environment e. We are now in a position to calculate 
the interference for the propagation ([2]). To obtain explicit results, we consider particular 
initial states for the environment. We start with a single spin in thermal equilibrium, and 
later generalize to several spins in thermal equilibrium. 



B. Interference in a quantum system coupled to a single spin in thermal equilib- 
rium 

Consider the situation where the environment is a single spin of size (d — l)/2, which 
corresponds to a Hilbert space of dimension m = d. We assume that the energy levels of the 
spins are equally spaced, with neighboring levels separated by an energy HQ, as is the case 
for atomic or nuclear spins under linear Zeeman effect in an external magnetic field. In its 
own eigenbasis, the matrix elements of the spin Hamiltonian reads 

H^ = mv8 vp (3) 

where 1 < v < d and t/ — 1 is the number of excitations of the spin. We choose the spin 
to be initially at thermal equilibrium at temperature T = such that its density matrix 
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can be written as 



with partition function 



e 



-0H l 



tr(e-^) up Z up 1 } 



Z^(*) = £e-^ = i-^, (5) 



and x = (3 HQ. The propagator P simplifies 

d 

Z 



Pg/3,yS — Uap,"/uUp^Su e (6) 



Inserting into ([Tj), we finally obtain the expression for the interference in the propagation 
of S alone, 

n _ n d 

1 = V IP J 2 = — V I V U U* x e~ xv \ 2 

a, 7^<5 a, 7^5 



— P - x ( v + a )Tj 77* 77* /7 r 

/ / / j c ^ a/i,7f u ap,7ff u 

a, 77^(5 fi,v,p,a 



We are now in the position to investigate the statistical properties of I based on the statistics 
of U. Without prior knowledge of a particular set of quantum algorithms or physical time 
evolution, it is natural to choose U uniformly distributed with respect to the Haar measure 
dU of the unitary group U(N). The statistical ensemble for the joint propagator of system 
and environment is then the well known "Circular Unitary Ensemble" (CUE). This allows 
us in particular to recover previously known results [15| for the interference statistics for 
unitary propagation of S in the limit where the dimension of the environment is reduced to 
one, as we will show below. 

C. Numerical results 

For small dimensions n and m, one can obtain the entire distribution of interference 
P(X) numerically. We have produced numerically unitary matrices of size N = nxm drawn 



from CUE using Hurwitz parametrization 
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231 ] . In order to obtain good statistics we 



have used 10 6 matrices for the calculation of the distribution. Figure [T] shows P(I) for 
systems with sizes from n = 2 to 4, coupled to an environment of size m = 1 to 4 at inverse 
temperature x — 0.1. In the case n = 2, where the analytical calculation is possible the 




FIG. 1: (Color online) Numerically calculated interference distributions on a log-log scale for 
x = 0.1 for n = 2, 3, 4 (from the left to the right). In each plot, m = 1, 2, 3, 4 from the right to the 
left (red, blue, green, purple, respectively). The number of realizations is n r = 10 6 in all cases. 




FIG. 2: (Color online) Fit of numerically calculated P(I) (red curves) to a log-normal distribution 
(black dashed curves) at x = 0.1 for (n,m)=(4,2), (4,4) and (8,2) (from left to right). All fits are 
for n r = 10 6 except for the first plot in which the blue curve is for n r = 10 7 . 



distributions are very wide (see [15[ for m = 1 where -P(X) = 2 ^— ^ ). For higher values of 
n, the distribution becomes more and more peaked, and, on a log-log scale, more and more 
symmetric with respect to the maximum. The tails of the distribution decay more rapidly 
in the non-unitary case. We see that both the most probable value of X and the width of the 
peak decrease when m increases. As expected, the decoherence due to the coupling to the 
environment destroys the interference more efficiently with increasing m. For fixed m the 
general distribution behaves qualitatively like the distribution for the unitary case studied 



in [151 ]. i.e. the most probable interference and mean interference increase with n, whereas 
the width of the distribution decreases with increasing n. A change of the temperature 
essentially shifts the distribution. This is due to a change of the average interference with 
the temperature as we will see later (see eq.flTTJ)), and justifies why we have plotted all 
distributions for x = 0.1. 



S 



Figj2] shows that P(T) is well fitted by a log-normal distribution, 



exp 



/ (log(T)-M) 2 
y 2a 



PC) = ^7= L • (7) 



The fits work particularly well close to the center of the distributions, whereas deviations 
appear in the wings of the distribution. In addition, the wings appear to be clipped, but 
this is at least partly an effect of the finite number of realizations available. This is visible 
from the example n = 4, m = 2, where we have increased the number of realizations n r from 
10 6 to 10 7 . In the latter case, the clipping appears at substantially larger values of X. 

The numerically obtained distributions suggest that -P(X) is for n > 2 well characterized 
by its first and second moments. We will now present analytical results for these two 
moments which confirm the qualitative observations above for arbitrary values of n, m, and 
x, and make them more quantitative. 



D. Analytical results 



1. Average interference 



The average interference (X) follows from eq.( 



(I) 



-x (y-\-a) 



(8) 



where ( .) = J dU(.) means average over CUE. For the monomials composed of a relatively 
small number of factors U a ^ nv to be averaged here, the technique of invariant integration is 



well suited. We use the diagrammatical language introduced in [2J, |25|] to express (X) as 



X 



d 

ot,~1^8 ij,,v,p,(J 

Ctfl 



-- 8a 



-y 




(9) 



We refer the reader to 
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25 



261 ] for a detailed explanation and derivation of this technique, 



but summarize here the main features. For the sake of clarity we revert momentarily to single 
roman indices i, j etc. for rows and columns. All distinct row (column) indices that appear 
in the matrix elements of the monomial are represented by vertices on the left (right) with 
the corresponding label, irrespectively of whether or not they arise from a matrix element 



Uij or its complex conjugate U*j. A complex conjugate factor U*j is then represented by a 
thin solid line between the vertices i and j, whereas a factor Um is represented by a dotted 
line between the vertices k and /. When a given matrix element occurs with multiplicity t, 
a single line is drawn with the number t next to it to keep track of the multiplicity. Factors 
like \Uij\ 2 are represented by thick solid lines, which can also have a multiplicity larger than 
one. In 24j it was shown that the invariance of the Haar measure under arbitrary unitary 



transformations leads to the following important properties: 

(a) The value of a diagram does not depend on the specific values of the vertices. It only 
depends on the form of the diagram. This means that diagrams can be drawn without 
specifying the explicit values of the vertices. For example, 




(UnU^UuU^) = (u 24 u; 4 u 26 u; 6 ) = <r . (10) 



(b) If for at least one vertex in the diagram, the number of thin solid lines that originates 
from the vertex differs from the number of dotted lines then the value of the diagram 
is zero. For example, 

(u n u* 2 u; 3 u M ) = = o (ii) 

(UnUyinU^) = • ,/0. (12) 

In eq.Q, we sum over all row and column indices and different type of diagrams therefore 
appear, depending on which vertices coincide. Combinations of indices contribute for which 
the vertices (ju) and (5u) collapse on the vertices (ja) and (Sa), respectively, i.e. configu- 
rations with v = cr. We thus have 

n d 



'/. - * — ' * — ' •■•/' ^ i v 



n d d Su d 



«/x ^ \, 7^ 



At this point only two types of diagrams remain, and since their values do not depend on 
the summation indices, we get 
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The values of the two diagrams are easily found 24| , 

1 




N(N + 1) 
1 

N(N 2 - 1) ' 



The prefactor can be rewritten as 



h(x) = Z }^ ( X \ = coth(rfx/2) tanh(x/2), (16) 
Z (x) 



and we finally obtain, with d = m, 

,mx. , ,i.nm(n-l) 2 
— tanh- ) \ ' 
2 2 (n z m^ — 1 



(J(n, m, a;)) = coth( )tanh(-) . Q v V- ( 17 ) 



This is our first central result which we now discuss in detail. 

We first observe that the entire temperature dependence is entirely contained in the 
prefactor h(x). Its limits for x — > and x — > oo are 1/m and 1, respectively. In between, h(x) 
increases monotonously. We thus find that the average interference decreases with increasing 
temperature, an intuitively appealing result. Only the dimension of the environment m = d 
enters the dependence on temperature. This is true in fact for all moments of -P(X), as the 
entire temperature dependence is contained in factors exp(— xv) which are always summed 
over v = 1, . . . , m. 

In the particular case m — 1, i.e. n = N, we recover as expected the expression for purely 
unitary propagation fis| . 

{I{n , hx)) = ^f = ^Jl = {Inm . (18) 

No entanglement or correlations with the environment can arise in this case, as a single state 
always factors out, such that the dynamics of S remains indeed entirely unitary. 

Contrary to what might be expected naively, the unitary result is not recovered for zero 
temperature, x — > oo. Rather one finds 

j\[( n - l) 2 

lim(X(n,m,x)> = \ ! , (19) 
We recall that N = n x m. For n ^> 1 and m fixed we have the asymptotic behavior 

(Iu(n)) = n -2 + 0(-) (20) 

n 

(l(n, m, x - oo)) = ^— - + O(-) ~ (21) 

m n m 
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We see that for n 3> 1, the average interference still scales linearly with the system size, 
but is roughly a factor m smaller than in the unitary case. The reason for this reduction 
is, of course, that even for a heat bath initially in a pure ground state, the common uni- 
tary dynamics of S and S entangles S and S, such that after tracing out the environment 
non-unitary evolution of S results. The consequent loss of coherence manifests itself in a re- 
duction of interference. In the opposite limit of infinite temperature, x — > 0, the temperature 
dependence of the prefactor h(x) leads to reduction by another factor m, 

hm(J(n,m,x)> = ^— j . (22) 

The additional reduction is also seen in the asymptotic expansion for n>l, which reads in 
this case 

In, m,x = = — =- + O - ~ V U \ J/ . 23 

For m ^> 1 and n fixed we find 

(J(n, m, x - oo)> = (n ~ 1)2 + 0( Ar) (24) 

urn rrr 

(n-1) 2 1 

(J(n, m, x = 0)) = K - J- + O(-) . (25) 

Eqs. (I24j) and ( 1251) show that for fixed n > 1, (X) decreases as 1/m (1/m 2 ) for zero tempera- 
ture (infinite temperature). In Figj3]we plot (X(n,m,x)) for four different temperatures as 
function of n and m. We see that for given temperature, (X(n,m,x)) increases with n, but 
decreases with m. For large n, with m and x fixed, the increase is essentially proportional to 
n, just as in the unitary case, albeit with a slope reduced by a factor h(x)/m. For large m, 
with n and x fixed such that mi 3> 1, the decrease of (I(n,m,x)) is roughly as 1/m with 
a prefactor "^"^ ■ More generally, an increase in the dimension of the environment 

decreases the average interference in a power law fashion with a power that crosses over 
from m~ 2 for x = 0.001 to m" 1 for x=10 and fixed n. One should not conclude from this, 
however, that a quantum system coupled to an infinite dimensional heat bath will never 
show any quantum interference effect. Rather, it should be kept in mind that we consider 
here generically strong couplings to the environment, in the sense that a typical joint evo- 
lution operator U of S and £ does not distinguish the two subsystems, or, for that matter, 
a system hamiltonian, bath hamiltonian, and coupling hamiltonian. It is natural that such 
strong couplings destroy coherence and thus quantum interference rapidly, but the situation 
can of course be different for weak couplings. 

12 
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FIG. 3: (Color online) Contour plot of ln((T(n, m))) for x=0.001, 0.01, 0.1 and 10 (upper left to 
lower right), for n and m between 2 and 1024. The distance between the contours is 2, and the 
red dashed line is for ln((X(n, m))) = 0. Values increase from dark to bright colors. 
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2. Second moment of the interference distribution 

In order to appreciate the width of the interference distribution as function of m, n and 
x, we now calculate the second moment of -P(X). By taking the square of the eq.([7|) we find 

1 n n d d 

T 2 - — \ ^ \ ^ \ ^ \ ^ -x(u+U+n+S)jJ J J* J J* TJ TJ* TJ J J TJ* 

' ^4 / / / , / j / j e ° ap^v^ a p,6v u ap,~fo u ap,&(r u ap,gq u ap,dq^ ar,gs u a r,ds 

a, 7^(5 a,gj^d p,p,p,r v,a,q,s 

The fact that 8 factors U appear now, makes the analytical calculation of (I 2 ) rather cum- 
bersome. As we will see, altogether 19 different diagrams contribute. We give here a rough 
outline of the derivation, relegating most details and in particular the values of all dia- 
grams to the Appendix. In order to streamline the presentation we introduce the following 
simplifications of notation: 

• First, for the six subsystem indices, we substitute (a, 7, S, a, g, d) — > 

(«i, a 2 , a 3 , a 4 , a 5 , a 6 ). 

• Similarly, for the eight environment indices, we replace (//, p,p,r,u,a,q,s) — > 

(/^l, ^2, /^3; /^4, ^5) /^6) ^7) fl>8)- 

• We then drop the redundant letters a and /i altogether, both from matrix elements 
and the diagrams. I.e. we write matrix elements C/ aiWiQ!fcW just as Uij^i- So now Lq^n 
is not the first element of the matrix but it is the element with indices (ai/ii, 
Recall that all a (/i) indices take values between 1 and n (m), respectively. 

• The constraints 7^5 and g ^ d read now a 2 7^ a 3 and 7^ a§. They are assumed 
implicitly. 

• We also make it a rule that in a sum fij} denotes the set of all indices which 
appear explicitly in the summand as indices of matrix elements, or, equivalently, as 
labels of vertices, with the exception of those which appear under another sum in the 
same expression. E.g. in J2i a u \ a > the ^ YS ^ sum * s over an a ' s an d ai l A*' s that 
show up in the diagram summed over, with the exception of /is and /i7, which are 
considered separately. 
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We can then write 



(n,m) 

( J2 > = «E e^^^ + ^^)(^ n ,35^r i ,45^r 2 ,3 6 ^12,46^3,57^3,67f/24,58f/24,68) (26) 



z 4 



5.8 
6.7 



-2 (At 6 +M6+M7+At8) 2 ' 3 " ~ 5,7 (27) 

Z 4 ^ - 4.6 ' 




We re-emphasize that the indices of £7 which appear in eq. (l26i) are indices of indices, 
e.g. C/ii.,35 = U aifllja3fl5 . As for equation ([8]), the only non-vanishing contributions arise 
from diagrams without open ends. They correspond to three distinct configurations of the 
summation indices, namely /i 5 = /i 6 and /17 = /i 8 , or a 3 = a 5 , a 4 = a 6 , /i 5 = /i7, and /i 6 = /ig, 
or «3 = «6,a4 = a5,/i5 = /is, and /i6 = /i7- These three configurations give rise to three 
sums, 

, , 2,4 -^r — 6,7 2,4 -^r — 4,6 

(n,m) m (n,m) n m X 

<*> - u e E^ h+ " eee 

{o ij( Uj}M5,M7 "X" o c {aj,/^} 03,04^5^*6 EXE oc 

1,1 3,5 1,1 3,5 

2,3 — ^- 4,6 

(n,m) n m X' 

+ EEE«-* t " , :r;:)' <*> 

1,1 X 3,5 

The last two terms are equal as can be seen by exchanging the summation indices <-> a^. 
We are therefore left with 



, > 2,4 x 6,7 2,4 ■ 4,6 

(n,m) m X (n,m) n m X 

= J_( Y^ e -2x(/,5+M7) 2 ' 3 - - 5 - 7 +2 e -2x( W +/ia)l,2 -" "- 3,6 

z 4 V ^ ^ 1,2 - 4,5 ^ ^ 2,3 • 4,5 

1,1 -X„rX 3,5 1,1 XT XX 3,5 

= ij(il + 2fl). (29) 
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The terms A and B depend on 19 different diagrams which we calculate again by invariant 
integration. For A we have 

, , 2,4 ^ 6,7 

(n,m) m ;;<; 

A= y Y< * "J - (30) 
Z — ' ^— ' 1,2 ^----ps- 4,5 
K, W }M5,M7 ,.; 

2,4 — 6,5 2,4 — — 6,7 

(n,m) m ><' m 

y ( e -4x M5 2 ' 3 — ^ 5 > 5 + e -2^^ 5 + W )2,3 — ^-s,7\ , 31 > 

V 1.2 1.5 -^-^ 1.2 1.5/ V ' 

1,1 *.<:...>.>,. 3,5 1,1 3,5 



2,4 — r^- 6 2,4 

2,3 " ■ 5 / ^ ^ a™ 2,3 * ■ 5 



(n,m) 

E {?, Zir !< E e " 4 *" 5 ) + "tZT!(E e-"" + "")) 02) 



. 1,2 ■ 1 ^ 1,2 ^-—-^ 4 

1,1 ^C . S^ 3 1,1 3 



2,4 — 6 2,4 

2,3 5 / x 2,3 ^- 5 



(n,m) 

= E < 33 > 

1.1 i?<> 3 1,1 J^L 3 

The single indices on the right hand sides of the diagrams in ([3"2j) now decode only en's; 
/ii . . . /i4 still appear explicitly on the left column of vertices, \i§ and /is were chosen identical 
to yU 5 and fi?, respectively, and the latter two indices are still summed over. A "bar" 
vertex is a vertex which cannot collapse with a "normal" vertex even if both values of the 
corresponding a's (or /x's) are the same. The vertices labeled 5 and 6 (which stand here 
for 0:5 and olq) inherit this property from the \x-j part still present in ([3"Tj) : The restriction 
//g /i7 implies indeed that none of the two top vertices can collapse with either of the two 
bottom vertices in the second diagram. Thus "normal" and "bar" vertices can only collapse 
on vertices of the same kind. The constraint ^ a§ is still implicit. The functions f(x) et 
g(x) are defined as 

m = (^3rJ = z{Ax) (34) 

fl_ e -2xd\ /i_ -2x(d-l)\ 

9(x) = e"- (y^r) ( j _ ^ ) = ~ ^) . (35) 

With this, we have introduced all notational innovations which allow the analytical cal- 
culation of (I 2 )- The rest of the calculations amounts to identifying all possible non-zero 
configurations of collapsing vertices allowed by the remaining summation variables. The 
explicit expansion of the terms A and B finally leads to 

(X 2 ) ^^[f{x)(A 1 + {n-l)A 3 ^+g{x)(^A 2 + {n-l)A A +n{n-l)B 1 +n{n-l) 2 B 2 ^. (36) 
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Here, the parameter d in the functions / and g is d = m. The terms A, and Bi are defined 
and calculated explicitly in the Appendix. They only depend on n and m. 




1234567 1234567 
ln(n) ln(n) 



FIG. 4: (Color online) Contour plot of ln(oj(n, m)) for x=0.001, 0.01, 0.1 and 10 (upper left to 
lower right), for n and m between 2 and 1024. The distance between the contours is 2 and the 
dashed line is for ln((X(n, m))) = —10, except for the last plot where the distance is 1 and the 
dashed line is for ln(aj(n,m)) = —6. Values increase from dark to bright colors. 

In FigJUwe plot the standard deviation of the distribution of X, ax(n,m,x) = ((X 2 ) — 
(X) 2 ) 1 / 2 for four different temperatures as function of n and m. For given temperature, 
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aj(n,m,x) decreases with n and m. The log-log-log plot shows that the decay behaves as 
a power law in n and in m. The corresponding powers can be found from an asymptotic 
expansion of the variance var(n, m, x) = a\ (n, m, x) for n ^> 1 or for m > 1 in the limits of 
zero or infinite temperature. For fixed m, we find for n ^> 1 

2(m-l) 2 4(m 4 - 3m 3 + 3m 2 - 5m + 3) n/ L 
var (n,m,x^ oo = ^ ^ ^ — '- + O — 37 

2(m 2 -l) 8 -4m 4 1 . 

mr n,m,i = = - ^ + =-=- + 0{— ) . 38 

This should be compared to the unitary case, where the asymptotic expansion reads 
varjj{n) = \ + 0(-\), as is still evident from eqs. (|3TI) and (1551) by choosing m = 1. We 
see that the variance decays more slowly as function of n in the presence of decoherence, 
i.e. as 1/n instead of as 1/n 2 in the unitary case. In other words, decoherence tends to slow 
down convergence of the interference distribution to a narrow peak. Nevertheless, the power 
law decay of the variance as function of n implies that, also in the non-unitary case, the 
interference distribution becomes for u>la very narrow peak centered about the average 
value (which itself increases with n, see eqs.( H9p and (T221 ). 

Asymptotic expansion of var(n, m, x) as function of m ^> 1 with fixed n gives 

2(n- l) 2 1 , 

v ar(n, m, x — ► oo) = = — 5 h 0( — 7 ) (39) 

(n-1) 2 1 

varin, m, x = 0) = ^ — j- + 0(— ) . (40) 

Thus, also an increase of the dimension of the environmental Hilbert space narrows the in- 
terference distribution. However, since according to (I24ll25p . the average interference decays 
as 1/m (1/m 2 ) for x — » oo (x — > 0), the relative width, i.e. standard deviation divided by 
the average value, is asymptotically independent of the dimension of the environment. 

In the case m = d = 1 (n = N) all the pref actors m[i] (see Appendix) are zero if 
i > 1. With the same parameters we have furthermore from eqs. (I34ll35j) . f(x) — Z — 1, and 
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g{x) = 0. Thus the expression (136|) of (Z 2 ) simplifies considerably. 
(J 2 ) = (NA 1 + N(N - 1)A 3 ) 
V(iV[3]^^ + 4N[2] 



2JV[1] 




) 



+iV(iV - 1) (JV[3] < ^ + 4A^[2] ^> + 2N[1] X ) 
N (N 3 - + 8) - 4 



(N + l)(N + 3) ' 

As expected this leads to the standard deviation erj 
for purely unitary propagation [H 



AT-1 



JV+1 V Af+3' 



identical to the expression 



The numerical results in section III CI are in very good agreement with our analytical 
results, as can be seen in table [J where we compare the numerically obtained average values 
and standard deviations for the examples shown in figfSJand for (n, m) = (4, 8) and (8,4) to 
the corresponding analytical results. 
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m 


{I) (num.) 


{I) (ana.) 


ox (num.) 


ax (ana.) 


4 


2 


0.57279 


0.57286 


0.11728 


0.11719 


4 


4 


0.14296 


0.14293 


0.03260 


0.03255 


4 


8 


0.03702 


0.03702 


0.00864 


0.00864 


8 


2 


1.54120 


1.54109 


0.09022 


0.09409 


8 


4 


0.38796 


0.38796 


0.02670 


0.02666 



TABLE I: Comparison of numerical and analytical values of (I) and ax- All results are rounded 
to five digits after the decimal point. 



III. INTERFERENCE FOR A SPIN COUPLED TO SEVERAL SPINS 

In this part, we generalize the previous calculations to a situation where the environment 
consists of s independent spins with d energy levels with energy spacing Kl. Thus, the 
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dimension of the environment is m = d s . The hamiltonian of this system reads 

ff« = £X>, (41) 

k=l 

where is the hamiltonian of spin number k (eq.Q). The components of in its 
eigenbasis are 

with the notation for the indices v — (ui, v 2 , v s ) and p = (p%, p 2 , p s ) 

The density matrix corresponding to the thermal state of such a system factorizes, e^ s ' = 
= e (!)® s ; which leads to the components 



-xS{v)s 



with x = (3hfl, S(v) = X/L=i u k, and where Z is the partition function of the thermal state 
of a single spin introduced in the previous section. It turns out that in order to generalize 
the previous calculation of < X > and < X 2 > to this kind of environment, we just have to 
replace Z by Z s in eqs. ( |T5l) . ( |34|) . and ( 1351) . and keep d = m s instead of d = m. This is again 
a consequence of the fact that the values of the diagrams do not depend on the indices of the 
vertices. Thus, the same values are obtained even for composite indices reflecting several 
subsystems, and only the multiplicities and temperature dependent factors are modified. 
Since the spins of the heat bath are taken as non-interacting, the sums over the thermal 
factors just gives rise to powers of the single spin thermal factors, as is the case also for the 
calculation of the partition function for s spins. This means that we have to replace 

d d 

f(x) -> Y, e ~ ixSM = Y,-Y. e ~ 4x ^- e ~ 4xfls = zs ^ = f s w 

M Ml Ms 

g [x) -> ^ e -M%)+sM) = ^ e -MSM+s(»)) _ ^ e -^W = Z s (2x) - Z s (Ax) . 
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The expressions for (I) and (T 2 ) become 



(J 2 ) 




(44) 



(45) 



+g s (x) [A 2 + {n- l)A A + n(n - l)5i + n(n - l) 2 B 2 



) 



(46) 



The argument m in the terms A, and i?j in eqs. (144(4*61) is now m = d s . It means that s 
spins of size (d — l)/2 act very similarly as a single spin of size (d s — l)/2, when it comes 
to their influence on the first and second moments of -P(X). The only difference lies in the 
temperature dependent prefactors f(x),g(x) and h(x). For a single spin of size (d s — l)/2, d 



have s = 1 for a single spin. For s spins of size (<i — l)/2 the dimension d in eqs. (U6ll34|[35|) 
remains, and s is the number of spins in eqs. (I44II46P - In the limits x — > or x — » oo the 
expressions coincide for the two situations. 

IV. SUMMARY 

We have investigated quantitatively how quantum interference is affected by decoher- 
ence. Based on a distribution of unitary matrices drawn from CUE which describe the joint 
propagation of system and heat bath, we have shown that the average interference increases 
roughly linearly with the Hilbert space dimension n of the system, but decays as a power 
of the dimension m of the environment. That power depends on the temperature of the 
environment (chosen here as one or several non- interacting spins), with a decay that essen- 
tially scales like 1/m 2 for T = 0, and as 1/m 3 for T — > oo. The width of the distribution 
decreases more slowly when decoherence becomes important, but for fixed m, the width of 
the distribution still decays as l/\/n (instead of as 1/n in the unitary case). Thus, for n 3> 1 
and m fixed, the distribution of quantum interference is still a sharp peak concentrated on 
the average value. Numerically we have shown that the interference distribution in the non- 
unitary case can be well fitted to a log-normal distribution for sufficiently large n, which 
implies that the number of i-bits [3] is to good approximation Gaussian distributed. 



m eqs. 



is given by the dimension of the environment m = d s , but in eqs. (|44f46l) we 
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V. APPENDIX 

We provide here the remaining details of the calculation of the terms A and B in the 
expression for (X 2 ), eq. (|29|) . as well as the values of the resulting diagrams. 

A. The A term 

^From eq. ( |33i) we have 



(n,m) 



.4 



1,4 
1,3 



E ( E mZZZ^'K 

{ai,fij} { ai =a 2 } 



1,1 



1,4 
1,3 
2 
1,1 



{0:1^02} 



2,4 
2,3 
2 
1.1 



2,4 — 
2,3 



4 ^ V ' 1,2 



(n,m) 



;r 



i.i 

4 - 

>r 

1 



+n{n — 1) /(x) 



4 -=rr — — 6 4 

^ " + ™0 - 1) 9(x) I 

2 4 2 

1 3 1 

= n + ng(x)A2 + n(n — 1) f(x)A^ + ^(n 

By taking into account the constraints on the aij we get 

4 — 6 



l)g(x)A 4 



(n,m) 

E 



E 



/ 4 

3 



n[3] 



V 



+ 4n[2] 



+ 2n[l] 




n[3]A 11 +4n[2]A 12 + 2n[l]A 13 . 



(47) 
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withn[i] — n(n—l)(n—2)...(n—i). We check that we have the n[3]+4n[2]+2n[l] = n 2 (n-l) 2 
configurations corresponding to the sum over the four indices aj with the two constrains 
a 3 7^ a 4 and a 5 ^ a 6 . The A lk read 



in 



E 



m[3] • m 2 (Y -2 ) 



+m[l] (V + + 4^g^) +m 




112 



E 



m[3]^?^ +2m[2] 




+ 




+ 



l 13 



E 




+ + J + 4m[l] + m 






m[3] --"^ + m[2](4<^ - 2 : ) 

+m[i](2X + : . • 1 ^) + 

For A 2 we obtain directly 



m 



4 — 



(n,m) 

^ = E 2 




= n 2 (n-l) 2 A n , 
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whereas A% is given by 



(n,m) 

E 

{ai,Hj} 



4 
3 
2 
1 

/ 4 



= E 



n[3] 



V 



+ 4n[2] 



3 



n[3]A 3 i + 4n[2]A 32 + 2n[l]A 33 . 



The A 3fe are 



+ 2n[l] 




/ 



A 



31 



E 



= (to[3] + 4to[2] + 2m[l] ) + (2m[2] + 4m[l] ) + (to[1] + to) 



A 32 = ^ 



4 

3 
2 
1 



(m[3] + 4m[2] + 2to[1]) ' + (2to[2] + 4m[l]) + (m[l] + to) ^> 



^33 = E 

{Mj} 





= (ro[3] + 4to[2] + 2m[l] ) + (2 m [2] + 4m[l] ) <t§C + ("' 1 + "') X 




The term A 4 can be expressed in terms of A 31 , 

4 jt- 6 



(n,m) 

■44 = E 



n 2 (n- I) 2 A 



31 



As a consistency check, we verify in the calculation of the terms A 3i that we have the 
to [3] +6to[2] + 7to[1] +to = to 4 configurations corresponding to the sum over the four indices 
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B. The B term 



In the same way as for A, we find for the B term 

, , 2,4 - 

(n,m) n m 

* = E E 5> 



-2x (/i B +/ie) 1,2 



2,3 
1,1 



(n,m) 

= E E (**) 

{oi,/i-,} 03,04 



2,4 ^- 4 

1,2 
2,3 

1,1 • 3 



(n,m) 

n(n-l)^(x) ^ 

{Oi,Mj} 



2,4 
1,2 
2,3 
1,1 



= n n 



it 1 ]} 01 02 



2,4 - 
1,2 
2,3 
1,1 



= ram 



ra 



+ n(n 



= ra 2 (ra- l)p(rr) (b y + ( n - 1)S 2 ) , 
where the terms Bi are given by 



4 — 
3 

2 •< 
1 



111 



= m[3] + 4m [2] 



+ 2m [2] 



+m[l] 



+ 2J§C +4 




m 




m[3] +4m[2] + 2m [1] 



+ (2m [2] +4m[l] 
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C. Analytical expressions for all diagrams 



All diagrams can be calculated by invariant integration. 
i -i 



We find 




jV(JV+l) 



(N+3)(N+2)(N+1)N 
N 2 +N+2 



(N+3)(N+2)(N 2 - 


-1)N 2 


-4 




(N+3)(N+2)(N 2 - 
-2 


■1)N 


{N+3){N+2)(N 2 - 
1 


-1)N 


(N+3)(N-1)N 2 
2 




(N+3)(N+2)(N 2 - 


-1)N 2 


4 




(N+3)(N+2)(N 2 - 


-1)N 2 


-(N 2 +l) 




(N+3)(N 2 -4)(N 2 
1 


-1)N 2 


(N+3)(N+2)(N 2 - 
2 


■1)N 2 


(N+3)(N+2)(N 2 - 


■1)N 2 


(N 2 +6) 





N(N 2 -l) 




1 





(Ar+3)(A r +2)(Af+l)Af 



(N+3)(N+2)(N 2 -1)N 2 

= N+l 

(N+3)(N+2)N 2 (N-1) 



^1 

{N+3)(N+2)(N+1)N 2 



^1 

(Af+3)(Af+2)(A f2 -l)A r 



(N+3)(N+2)(N+1)N 2 



3jV-l 

(N+3)(N 2 ~4)(N 2 -1)N 2 



(N+3){N+2)(N 2 -1)N 2 



~(N 2 +2N+2) 
(N+3)(N 2 -4)(N 2 -1)N 2 



(N+3)(N+2)(N 2 -1)N 2 





(N 2 -9)(N 2 ~4)(N 2 ~1)N 2 
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